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1. INTRODUCTION 

Recently, speed-over-ground (SoG) radar has attracted attention as a technique to estimate vehicle 
speed through contactless measurements [1]-[4]. Unlike the conventional wheel speed sensor, SoG radar can 
deliver true ground speed because its estimate is independent of wheel conditions such as slipping and 
locking. It is expected that SoG radar will become more prevalent in the future, owing to the increase in 
demand for true ground speed measurement for vehicle safety systems such as anti-braking systems (ABS) 
and electronic stability control (ESC) [5]. 

SoG radar technique utilizes Doppler frequency to estimate the vehicle speed. There are several 
studies based on spectral estimations [6], [7], with the most popular method is based on the fast Fourier 
transform (FFT) because it has low computational complexity and is easy to use. This method is based on the 
decomposing sinusoidal component of the Doppler signal whose phase and amplitude vary across 
frequencies. Several approaches can be used to compute the FFT power spectrum, some of the methods that 
have been evaluated for SoG radar are periodogram [8] and Welch’s method [9]. Nevertheless, the 
FFT-based method has a critical issue with SoG radar. In practice, the Doppler frequency is non-stationary 
and can change rapidly according to vehicle motions, but FFT requires the frequency within the analysis 
window to be stationary. To adapt to the rapid change of Doppler frequency, the analysis window is fixed to 
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a short length, with a typical duration of fewer than 100 milliseconds. However, such a short length produces 
poor spectral resolution which consequently affects the accuracy of SoG speed measurement. 

Spectral estimation based on the autoregressive (AR) model has also been evaluated for SoG radar. 
The choice of using the AR method is largely motivated by improved frequency resolution for short data 
length [10]. Unlike the FFT approach, this method models the data outside the processing length, which 
improves the spectral resolution of the observed signal [11]. AR spectrum can be computed using several 
techniques including the Yule-Walker and Burg algorithm [12]. However, the main disadvantage of the AR 
method is that they are more computationally intensive compared to FFT and incorrect modeling of the signal 
may lead to statistical instability of spectra. 

This paper investigates the performance of SoG radar designed based on two spectral estimators 
namely periodogram which is based on Fourier transformation and Burg algorithm which is based on AR 
method. For this purpose, the theoretical aspects of both Fourier transformation and AR were reviewed. A 
simulation approach is used to test the spectral estimator and to provide comparative data of performance. 
Finally, graphical plots and tables are used to explain the performance of each estimator. The remainder of 
this article is arranged in the following manner. The signal processing formulation and methodology are 
presented in section 2. The simulation procedure is presented in section 3. The results and their performance 
analysis are presented in section 4 and finally, section 5 concludes the work on this paper. 

The assessment of spectral estimator performance for SoG radar requires an understanding of the 
processes underlying the generation of Doppler signal spectrum, referred here as “Doppler spectrum’. 
Consider a SoG setup illustrated on the top part of Figure 1. A radar with an antenna beamwidth of @ is 
mounted on a vehicle and faced the ground surface at an angle of a to the direction of the moving vehicle 
with a velocity of v. The propagation of waves emitted from the radar to the moving ground surface 
eventually backscattered toward the radar. Accordingly, Doppler shift from each scatterer on the ground is 
accumulated in the radar as the vehicle moves longitudinally above the ground. The resultant spectrum has 
peak magnitude when the antenna boresight beams directly to the moving ground and the magnitude reduces 
as it goes away from the boresight. 


Power (dBm/Hz) 


i 


Frequency (Hz) 


Figure 1. The principle of Doppler spectrum generation 


The Doppler spectrum of the first half-power arising from that process has a shape that can be 
approximate by a Gaussian function [7]: 
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where G(f) is the spectrum coefficient, o is the standard deviation, and f and fo is the Doppler frequency, and 
mean Doppler frequency of the distribution. The standard deviation and the mean frequency of the spectrum 
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are related to the radar antenna beamwidth and the vehicle speed respectively [13]. The (2) and (3) provide 
the main characteristics shape of the Doppler spectrum. 


foa 
~ 2cosa (2) 
v . 
o~-Osina (3) 


The location of the mean Doppler frequency is at the peak magnitude of the spectrum. Nonetheless, 
the actual spectrum found in practice has random amplitudes near the peak of the spectrum [13]. This 
amplitude fluctuation is due to the random scattering processing of microwaves on the ground surface which 
randomly modulates the amplitudes of the spectrum [14]. Thus, making detection of mean frequency 
difficult. 


2. RESEARCH METHOD 

The proposed framework for evaluating the spectral estimators is shown in Figure 2. In this 
approach, the spectral estimators estimate Doppler spectra provided. Afterward, the mean Doppler 
frequencies are extracted. Finally, the spectral estimator performances namely bias, standard deviation, and 
root mean square error (RMSE) are recorded. All numerical simulations are performed in matrix 
laboratory (MATLAB). 


RMSE, bias and 
standard deviation 
estimation 


Doppler signal Power spectral Doppler mean 


dataset estimator frequency extraction 


Figure 2. The framework of spectral estimator evaluation 


2.1. Synthetic Doppler signal and dataset 

The simulated signal in this work is generated in the same way as in [15]. Its parameters are: vehicle 
speed, v mean Doppler frequency, fm radar antenna beamwidth, 0 radar beam depression angle, a carrier 
frequency, fe sampling frequency, fs and signal-to-noise ratio, SNR. Additionally, the yielded “bell-shape” 
spectrum such as shown in Figure 1 is incorporated into the model using (2) and (3). All Doppler signal in 
this study is created using the following radar parameters: a=45 degrees, 0=24 degrees, f-=24.125 GHz, and 
fs=25 kHz. The antenna radiation pattern is assumed to have a Gaussian pattern at the first half-power. 

Doppler signals were generated with a mean Doppler frequency range between 0.1 kHz and 4.0 kHz 
at 0.1 kHz intervals. For each mean frequency, 500 Doppler signals were generated using the same mean 
frequency but with randomized white noise. On top of that, Doppler signals were also generated for four 
levels of SNR values of 20 dB, 30 dB, 40 dB and 50 dB. A total of 80,000 Doppler signals were generated 
for the evaluation process. It should be pointed out here that the four levels of SNR value (20 to 50 dB) and 
range of Doppler frequencies selected in this work are based on some states of practical SoG operation found 
in practice. Table | indicates the relationship of SNR, and mean frequencies to some states found in the real 
world. Figure 3 shows an example of 500 samples of time-series Doppler signal for vehicle speed moving at 
constant speed v=40 km/h which corresponds to a mean Doppler frequency, fm=1 kHz, and SNR values of 
SNR=40 dB. The in-phase and quadrature components of the signal are indicated with a solid and dashed 
line, respectively. 


Table 1. The relationship between test values and actual state conditions 


Parameters Values State conditions [16] 
SNR 20-50 dB 20 dB: Water surface in excited condition 
40-50 dB: Dry asphalt surface 
Mean Doppler frequency, fm 0.1 to 4.0 kHz Vehicle speed between 2.9 to 126.7 km/h 
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Figure 3. An example of time-series Doppler signal 


2.1. Power spectral estimator 
2.1.1. Periodogram 
The periodogram of an N points of discrete Doppler signal, x(0), x(1), ... x(n-1) is given by (4) [17]: 


Pre(f) =+ EN x, exp( — j2nfn)? (4) 


where f is the frequency. The sequence of x(n) is divided into segments of fixed length window, and at the 
same time, a window function is applied to each segment before the periodogram computation. Typically, 
fast Fourier transform (FFT) is used for an efficient computation of DFT, in which the signal sequence is 
divided into frames of 2” i.e., 256, 518, and 1024. The equivalent frequency resolution is 1/N, for example, 
the periodogram sample sequence is 0, 1/N, 2/N...(N-1)/N. Since the frequency resolution is inversely 
proportional to the sample length, the precision of frequency measurement is restricted by the frequency 
resolution. A windowed sinusoid in the time domain produces sin(x)/x function with sidelobes in the 
transform domain. Depending on the windowing function used, the sidelobes level can have a substantial 
effect on the effective Doppler signal SNR of the transform. 

The number of periodogram frequency samples can be increased by the interpolation process. 
However, this approach is computationally costly. An efficient approach is to perform zero-padding; a 
process of which adding zeroes in the discrete Doppler signal before the FFT computation. 


x(n)ifn < N 
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where M is the length of the zero-padded signal and must satisfy the M>N requirement. It is worth 
mentioning that the zero-padding process does not add more information nor increase the spectral resolution 
but, artificial increases of frequency samples, thus make it appear that the frequency resolution is increased. 
Nonetheless, the interpolated periodogram samples may help to reduce the bias in the estimation of the peak 
frequency component. 


2.1.2. Autoregressive (AR) 
For a given Doppler signal x(n), the AR method estimates the model of the Doppler signal as an 
output of a linear system driven by white noise [17]. 


x(n) = — Leas p(k) x(n — k) + e(n) (6) 
Where a, is the AR coefficient, p is the AR model order, and e(n) is the white noise. The AR method 
produces a higher spectral resolution than the periodogram approach for short records because the PSD is 
computed from the model of the Doppler signal. Additionally, the AR method can also resolve two close 


sinusoids if the model is accurate, and the SNR is high. 
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The procedure of estimating power spectral density (PSD) based on AR consist of two main steps; 
first the AR parameters are estimated from the sequence of Doppler signal, x(n), and secondly, the PSD is 
computed. There are many methods to estimate the AR parameters such as Yule-Walker and Burg 
Algorithm. In this paper, we use the Burg algorithm because it produces a stable AR model, and is relatively 
efficient [18]. The PSD estimated using the Burg algorithm is defined by (7): 


Pap (f) = i (7) 


i Z 
1+2? ap(k)e-j2nfk]| 


where ep is the total of the least square error that computed from the sum of the forward prediction error fp 
and backward prediction error bp. 


fon) = x(n) + XÈ ap(k) x(n — 1) (8) 
b (n) = x(n — p) + Xk- ap (k) x(n + 1) (9) 


The details on how the Burg algorithm estimates the AR parameters can be found in [19]. The AR 
model order can be any order as required with the spectrum approximation become increasingly accurate 
with the increase of order [20], [21]. Nonetheless, the computational cost increases with the increase of 
model order with the number of computational operations is proportional to O(p°) [22] where O is the 
operation complexity. Furthermore, high model order increases the variance of the spectrum [23]. A suitable 
p-value that can sufficiently exhibit the properties of the signal but is not too large can be determined using 
the Akaike information criteria (AIC) and it is given by (10): 


AIC(p) = ® + log( 04°) (10) 
where g4” is the variance of the estimated prediction error of white noise of model order p. 


2.2. Mean frequency estimation 

We used two methods namely the maximum magnitude of the smoothed spectrum, and the spectral 
centroid method to estimate the mean frequency from the Doppler spectrum. They were different on the basis 
that they differ in terms of approach. Firstly, the maximum amplitude is based on selecting the maximum 
magnitude of a smooth spectrum. To reduce the effect of random magnitude fluctuations about the spectrum 
peak, the maximum magnitude is selected after the spectrum is average with an average moving filter of 
n=10. Secondly, the spectral centroid is based finding the mass center of the spectrum [24]. 


Centroid = Mass Lease (11) 
Where f(n) is the frequency bins and S(n) the power of the frequency bin. To reduce the effect of spectral 
noise into the centroid computation, the center-of-mass is determined by selecting the spectrum main lobe 
frequency components. In this work, the region of the spectrum main lobe of significant SNR is isolated from 
the rest of the spectral components. The isolation is performed by removing the components that are less than 
10 dB from the noise floor. 


2.2. Accuracy indicators 

A comparative assessment of the spectral estimators is made with the use of three statistical 
performance indicators. They are bias, standard deviation, and root mean square error (RMSE). 
Mathematically, they are calculated as (12)-(13) [25]. For every estimation of mean frequency from Doppler 
signal in the dataset, the accuracy of the spectral estimators are measured using a set of core statistical 
performance indicators namely bias, standard deviation, and root mean square error (RMSE), and they are 
defined as (12)-(13) [25]. 


Bias = f — fin (12) 


StdDev = |——YN=3(fn(n) — FY? (13) 
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RMSE = [EINAR = fr)? (14) 


Where fm the measured mean frequency, fn is the true value of mean frequency and f,, is the average mean 
frequency of N observations. 


3. RESULTS AND DISCUSSION 

Figure 4 visually indicates the crude accuracy of approximating mean Doppler frequencies as a 
spectrum peak using the periodogram and AR-Burg method with Figures 4(a)—(1) show different combination 
of mean frequency and SNR. The spectra for periodogram and AR are presented using solid and dashed lines 
respectively. The Ist, 2nd and 3rd row of the plot matrix shows the estimation for mean frequency of 
Jm=0.3 kHz, 2.1 kHz, and 4.0 kHz respectively with the true location is indicated with vertical dashed lines in 
the plots. The spectra were computed using 2048 samples of Doppler signal sliced using Hanning window. 
The equivalent period is 80 milliseconds, and the obtained frequency resolution is 12.2 Hz. For the AR 
method, the model order was computed for orders between 1 and 50 with order p=6 selected for all spectra 
presentation. As can be seen from the figure, despite the noisy appearance of the periodogram spectra, their 
peak is discernible even under the random amplitudes. This method is comparatively reliable to retain the 
shape of the spectrum main lobe across the value of Doppler frequencies and SNR levels which we consider 
crucial in the context of extracting the mean Doppler frequencies. 
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Figure 4. Power spectral density estimated using periodogram and AR Burg algorithm indicated with solid 
line and dashed line respectively, with (a) to (1) demonstrate various combinations of mean frequency and 
SNR 


On the other hand, for the AR method, the peaks were poorly approximated in some conditions. 
There were 2 cases in particular; in the first case, the AR method shows a single peak but somewhat shifted 
to the left from the true mean frequency for location fm=0.3 kHz and SNR levels between 20 and 40 dB. In 
the second case, the AR spectra appear to have split peaks for f,=4 kHz and SNR between 40 and 50 dB. 
These exhibits indicate the matching capabilities of the AR method to adapt with shape and SNR levels of the 
Doppler spectrum. 
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Note that for demonstration in Figure 4, the AR model order for all spectra is fixed at p=6. However, 
this value is not an optimal choice for all the Doppler signals in the dataset. We found that for AR method 
can approximate spectrum with a single peak for order range between p=3 and p=30 with peak become 
increasingly broader with higher model order. We have shown that for a certain combination of Doppler 
frequencies and SNR values, the AR spectral approximations were somewhat inaccurate. In these instances, 
the comparison of the spectral estimation methods is meaningless in terms of performance. Thus, for this 
reason, we considered some changes in the parameters of spectra would improve the spectrum estimation. 

In the case of fm=0.3 kHz, a simple approach to obtain an AR spectrum with a proper peak is to 
increase the model order. For example, a reasonably accurate single-peak formation can be obtained for 
fm=0.3 kHz and SNR=20 dB using model order p=15. However, this approach increases the computational 
cost as discussed in section 2 and therefore is not favorable for the low-cost design of SoG radar system. 
Instead of increasing the AR model order, we reduced the sampling rate. Our investigation found that a peak 
can be obtained when we reduced the sampling rate f=25 kHz by a factor of 5. However, this approach also 
reduces the time resolution to 0.41 seconds. For the case of split peaks found on fm=2.1 kHz and 4 kHz, we 
keep the same sampling rate, but we reduced the model order to p=3 to obtain a spectrum with single peak. 
Figures 5(a)—(1) show the AR spectra with different combination of mean frequency and SNR, for fm=0.3 kHz 
with new sampling rate, f=5 kHz with model order p=6, and AR spectra for fm=2.1 kHz and 4 kHz with 
standard sampling rate but with new model order p=3. 


SF, = 9.3 kHz, SNR = 20 dB Fi, = 9.3 kHz, SNR = 30 dB Í „7 9-3 kHz, SNR = 40 dB Fi, = 9.3 kHz, SNR = 50 dB 
0 0 0 
(b) (c) 
N N N N 
T 20 ME, E 20 T 20} fh z 
faa fes] fas # fes] 
© l kz) © l 4s) 
Pam AL F aan a 
z -40 i 2 -40 B -40 | A 
o io) io) e] 
a | a a l a 
| | 
-60 — -60 -60 -—— 0 sps 
0 500 1000 1500 2000 0 500 1000 1500 2000 0 500 1000 1500 2000 0 500 1000 1500 2000 
Frequency (Hz) Frequency (Hz) Frequency (Hz) Frequency (Hz) 
S,, = 2-1 kHz, SNR = 20 dB Soy = 2-1 kHz, SNR = 30 dB Ff, = 2-1 kHz, SNR = 40 dB Sin = 2-1 kHz, SNR = 50 dB 
0 0 7 
T Í N N N 
= .20 i = -20 | T T 
fes] | faa] | fes] fes] 
g B Zz B 
5 i 5 l 5 5 
> -40 i B -40 i z z 
o io) jo} io) 
a l a | Ay a 
l | 
-60 ; -60 l 0 -60 
0 1000 2000 3000 4000 0 1000 2000 3000 4000 0 1000 2000 3000 4000 0 1000 2000 3000 4000 
Frequency (Hz) Frequency (Hz) Frequency (Hz) Frequency (Hz) 
J „7 4.0 kHz, SNR = 20 dB SF, = 4.0 kHz, SNR = 30 dB Ff, = 4.0 kHz, SNR = 40 dB J „7 4.0 kHz, SNR = 50 dB 
0 0 
(i) 
N N N N 
= 20 - ~ T 40 i T T 
fes] fes] i faa [2s] 
gB 2 B 2 
5 5 i 5 5 
= -40 z -40 i 2 ž 
o e] io} e] 
a a l a, a 
l 
-60 -60 l -60 ! -60 ” . 
0 5000 10000 0 5000 10000 0 5000 10000 0 5000 10000 
Frequency (Hz) Frequency (Hz) Frequency (Hz) Frequency (Hz) 


Figure 5. Spectra of periodogram and improved PSD approximation for AR method with (a) to (1) 
demonstrate various combinations of mean frequency and SNR 


The average bias, standard deviation, and RMSE for 500 estimates of mean frequency for 
fm=0.3 kHz, 2.1 kHz, and 4 kHz as a function of SNR is shown in Table 2. For the periodogram method, the 
results were divided into two categories: estimation with maximum magnitude and estimation center-of-mass 
method. The recorded data are presented in relative percentage values to present the obvious level of 
accuracy for values of mean frequencies. Table 2 indicates that using a suitable method is important for 
certain performance criteria. For average bias, the periodogram method used with maximum magnitude 
produces the lowest bias. For variation in estimate and RMSE, the periodogram with center-of-mass and AR 
method produced the best result. As far as statistical accuracy criteria are concerned, it is difficult to conclude 
which estimators provide the best performance since each method has their strength. 
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Again, it is worth noting that for f,=0.3 kHz, The AR method was only successful if the spectra 
were computed with model orders p=15 or using lower frequency resolution. As for AR spectra for 
Jm=2.1 kHz and 4 kHz, the spectral estimates with a single peak were possible with p<4. This indicates that 
the AR method is guaranteed not to work for lower orders. Nonetheless, for a practical SoG radar system, 
this is not a major limitation since the many modern microprocessors are equipped with sufficient 
computational capabilities to support AR computation with high model order. The performance trends of 
each method as a function of Doppler frequency are shown in Figures 6(a) to (c). Although the estimated data 
fluctuate across the Doppler frequencies, we observed trends and we visualized the trend by using the linear 
regression model. The solid lines with markers correspond to the SNR values. 


Table 2. Mean frequency estimation using maximum magnitude. 
Periodogram and Maximum Magnitude Periodogram and Mass Centre Autoregressive and Maximum Magnitude 


A 


cae m Rel. Bias Rel. oo (%) Rel. RMSE Rel. Bias Rel. õọ Rel. RMSE Rel. Bias Rel. oo (%) Rel. RMSE 
kH) (%) di (%) (%) (%) (%) (%) Sis (%) 
50 0.3 0.3 7.1 7A 1.3 3.2 3.4 2.2 1.3 2.6 
2.1 0.2 6.1 6.1 1.1 2.5 2.7 2.8 1.4 3.1 
4.0 0.0 6.1 6.1 0.8 2.2 2:3 2.7 1.0 2.8 
40 0.3 0.4 7.1 7A 1.3 3.2 3.4 1.2 1.3 1.8 
2.1 0.2 6.1 6.1 1.1 23 2,7 1.3 1.3 1.8 
4.0 0.0 6.1 6.1 0.8 2.2 2.3 2.5 1.0 2771. 
30 0.3 0.8 6.9 6.9 1.2 2.8 3.0 1.1 1.3 1.7 
2.1 0.6 6.5 6.5 1.6 2.9 3.3 3.0 1.5 3.3 
4.0 0.5 6.7 6.7 0.4 2.0 el 1.1 1.0 1.5 
20 0.3 0.4 75 7.5 1.2 3.3 3.5 8.4 2.4 8.7 
2.1 0.4 5.5 5:3 1.4 2.3 2.7 8.3 1.9 8.5 
4.0 0.5 TA 7.2 0.6 2,2 2.3 1.0 1.1 1.4 

*— SNR = 50 dB —4— SNR = 40 dB —— SNR = 30 dB —+— SNR = 20 dB 


Ny 
© 


S 
© 


S 
© 


Periodogram + Max. Mag. Periodogram + Mass Centre AR + Max. Mag. 


n 
n 
n 


wn 
n 
wn 


Percentage Bias (%) 
S 
Percentage Bias (%) 
S 
Percentage Bias (%) 
S 


0 0 0 mi ws. 
0 1000 2000 3000 4000 0 1000 2000 3000 4000 0 1000 2000 3000 4000 
Doppler frequency (Hz) Doppler frequency (Hz) Doppler frequency (Hz) 
(a) 


N 
© 


Periodogram + Max. Mag. Periodogram + Mass Centre AR + Max. Mag. 


n 
n 
n 


Percentage Std. Dev. (%) 
5 

Percentage Std. Dev. (%) 
S 

Percentage Std. Dev. (%) 
S 


0 0 0 
0 1000 2000 3000 4000 0 1000 2000 3000 4000 0 1000 2000 3000 4000 
Doppler frequency (Hz) Doppler frequency (Hz) Doppler frequency (Hz) 
(b) 
= 20 => 20 => 20 
BS i Á r BS 
aed Periodogram + Max. Mag. = Periodogram + Mass Centre aad AR + Max. Mag. 
Ais E15 15 
n N n 
2 ž ž 
~~ 10 = 10 2 10 
50 Bp BD 
£ £ £ 
g5 35 35 
O o oO ó isiol 
5 3 Oo E lA ak, neces Dee 
a 0 ~ 0 a o : : 
0 1000 2000 3000 4000 0 1000 2000 3000 4000 0 1000 2000 3000 4000 
Doppler frequency (Hz) Doppler frequency (Hz) Doppler frequency (Hz) 


(c) 


Figure 6. The trend of relative percentage (a) bias, (b) standard deviation, and (c) RMSE of the spectral 
estimation methods 
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The estimated percentage of bias as a function of Doppler frequency is shown in Figure 6(a). The 
periodograms methods show the same linear trend which is apparent in both figures. There is a small increase 
in bias towards high Doppler frequency, but the increase is not appreciable. On the other hand, the AR 
method performance is somewhat non-linear, and they tend to have appreciable variation between SNR 
values. The AR method estimates were only relatively flat for Doppler frequency between the range of 1 kHz 
and 3 kHz but increased bias outside of the region, suggesting that model order may not be optimal. 
Subsequently, the estimated percentage of standard deviation as a function of Doppler frequency is shown in 
Figure 6(b). All estimates show an approximately linear trend, but between the two methods, the AR method 
shows the lowest increase in standard deviation with increasing Doppler frequency. Expectedly, the 
periodogram with maximum magnitude exhibited the largest variance in estimates. The estimated percentage 
of RMSE as a function of Doppler frequency is shown in Figure 6(c). Evidently, the increase in RMSE is 
accompanied by an increase in standard deviation. The periodogram method with center-of-mass method 
performs well indicated with the linear trend but with a small slope. As with the AR method, this method 
only shows a good prediction for mean Doppler frequency between 1 kHz and 3 kHz but deviated 
appreciably outside of this region. 


4. CONCLUSION 

This work compares two common spectral estimation methods namely Fourier transform based on 
periodogram, and AR method based on Burg algorithm to determine the mean frequency of Doppler signal 
found in SoG radar systems. A total of 80,000 artificial Doppler signals were generated and used as a test 
signal. The methods' performance was evaluated based on their measurement bias, standard deviation, and 
RMSE. The periodogram method is an effective method that can be applied to the time-domain Doppler 
signal even in the presence of considerable noise (20 dB). With the combination of mass-center detection 
method, the periodogram method can be said to provide good accuracy with consistent performance in the 
three criteria tested. Furthermore, this method is relatively simpler to use than the AR method. It can be 
deduced that, although the AR Burg method is theoretically superior in terms of resolution. It offered no 
significant advantages here. The AR method demonstrated only a slight advantage on RMSE when compared 
to the periodogram with center-of-mass method and was only appreciable within a narrow bandwidth of 
Doppler frequencies between 0.8 kHz and 1.3 kHz. Within this region, the AR method performs the best. 
Some factors such as the ratio of Doppler spectrum bandwidth to the sampling frequency might have a 
profound influence on the overall performance. The fact that fixing AR model order to a fixed value here 
seems to be a bad choice because they did not generate high consistency in the trends of the performance. 
This gives the impression that higher model order is desirable to better approximate the spectrum and thus, 
enhancing the performance of estimates. However, this study is not discussed here. 
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